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Abstract 

We have searched for the proton decay mode p — > v K + using the one-kiloton 
Soudan 2 high resolution calorimeter. Contained events obtained from a 3.56 kiloton- 
year fiducial exposure through June 1997 are examined for occurrence of a visible 
K + track which decays at rest into fi + u or 7r + 7r°. We found one candidate event 
consistent with background, yielding a limit, r/B{p — > v K + ) > 4.3 x 10 31 years 
at 90% CL with no background subtraction. 
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1 Introduction 



The prediction of many Grand Unified Theories, that nucleon decay occurs 
at accessible lifetimes, remains unverified but continues to motivate experi- 
mental searches. This expectation was foreshadowed in part by Sakharov's 
early suggestion that the simultaneous effects of baryon number violation, C 
and CP violation, and departure from thermal equilibrium could produce the 
baryon-antibaryon asymmetry observed in the universe. [jl] It is interesting and 
suggestive that no fundamental symmetry is known which implies the conser- 
vation of baryon number. Currently, nucleon decay as a consequence of the 
minimal SU(5) GUT model is considered to be ruled out experimentally[^],[|. 
However, other unification models, both with and without supersymmetry, 
predict baryon number violating processes. Amplitudes for these processes 
involve the exchange of new particles with unknown masses so that precise 
nucleon lifetimes are not predicted. The expectation that these masses will be 
in the range between the GUT scale of ~ 10 15 GeV and the Planck mass of 
~ 10 19 GeV leads to proton lifetimes in the range 10 32 — 10 35 years@]. De- 
cay modes with strange particles such as p — > v K + , are usually favored in 
models which incorporate supersymmetry |5|||. 

Previous searches for p — > v K + have been reported by the 1MB, Kamiokande 
and Frejus collaborations f?]fj,§]. The K + track can be imaged in ionization 
calorimeters such as Soudan 2 and Frejus, but is usually below Cherenkov 
threshold in water. 1MB searched for an excess of events in a region of anisotropy 
and energy with a large background |7[]. Kamiokande looked for an excess of 
single ring mu-like events between 215 and 255 MeV/c with a muon decay, and 
also for three-ring events compatible with an invisible, stopped K + — > 7r + 7r° 
decay[[|. Frejus used two-track events with ranges consistent with the K + and 
the /x+§. 

In the Soudan 2 analysis, we use both the visibility of the K + in a fine 
grained tracking calorimeter and the visibility of the decay electron from a 
stopped yU + to reduce backgrounds from atmospheric neutrino interactions. 
We searched for the proton decay mode p — > v K + using two K + decay 
channels, K + — > /U + z/ M and K + — > 7r + 7r°. 
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2 The Soudan 2 Detector 



The Soudan 2 detector is a time projection, modular iron tracking calorimeter 
with a total mass of 974 metric tons and fiducial mass of 770 tons. Details 
of module construction and performance may be found in References PHlll- 



The detector is assembled as a close-packed rectangular stack of 224 modules; 
each module is made from 1.6 mm thick sheets of corrugated steel, stacked 
in a hexagonal "honeycomb" structure. The average density is 1.58 g/cm 3 . 
On the walls of the underground cavern surrounding the detector, there is 
an active "veto" shield comprised of double-layer, hexagonal cell, aluminum 



proportional tubes[12 



Two million Hytrel plastic drift tubes (1.0 m long by 15 mm in diameter) 
fill the spaces in the honeycomb stacks. Ionization electrons deposited in an 
Ar/CC>2 gas mixture drift toward either end of the tube in a 180 volt/cm 
electric field with a velocity 0.8 cm//isec. Upon reaching the tube end, the 
electrons are detected by vertical anode wires and horizontal cathode strips. 
Each crossing of a tube by an ionizing particle can create an anode/cathode 
signal at a common drift time which we call a "hit". The pulse area, which 
is proportional to the integrated charge deposited in the tube, and the drift 
time for each hit are recorded by both the anode and cathode electronics. 
The primary trigger requires at least 7 hits, separated by at least 600 ns, 
in a group of 16 anode channels or at least 8 hits in a group of 16 cathode 
channels within a 50 /zsec window. The trigger efficiency for proton decay final 
states considered here is ~ 85%. The complete detector triggers at a rate of 
« 0.5 Hz from naturally occurring radioactivity and cosmic ray muons. Every 
240 seconds a "pulser" trigger provides a snapshot of the random background 
levels in the main detector. These are used as underlying events to add detector 
noise to Monte Carlo events. 



3 Event Analysis 



3. 1 Overview 



The data analysis proceeds in three stages. First we identify "contained events" . 
Event prongs are defined by scanning as track-like (n, fi or p) or shower-like 
(e or 7). Contained events are defined as having no hits on tracks or the main 
body of showers which are less than 20cm from the outside surface of the de- 
tector and the prongs do not start or end between modules. This is the same 
contained event selection as was used for our atmospheric neutrino analysis 
|0. Studies in reference |TB[ showed that the efficiency for correct identifica- 
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tion was 98% for tracks and 94% for showers. An absence of shield activity 
was required. Second, the events are required to have a topology consistent 
with the proton decay channel under study, based on counting the number 
of visible tracks and showers. Finally, kinematic selections which character- 
ize a particular proton decay mode are applied to the data and also to event 
samples which monitor background processes. 

The analysis procedure involves finding efficient selection criteria using our 
proton decay Monte Carlo program, while minimizing the backgrounds from 
atmospheric neutrinos and atmospheric muons. The former backgrounds are 
calculated using the atmospheric neutrino Monte Carlo program described in 
Reference [|13||, which incorporates the flux predictions of Barr, Gaisser and 



Stanev [|i4|| . Backgrounds from atmospheric muons may result when muons 
inelastically scatter in the rock outside the active shield. We use the term 
"rock" event to describe the interactions of a resulting secondary such as a 
neutron or Kl which goes into the Soudan 2 calorimeter and causes a contained 
vertex event. Most rock events have hits in the Soudan 2 shield which are time- 
coincident with the contained event in the calorimeter. These shield-tagged 
events are used to estimate any background from rock events without shield 
tags. A detailed analysis of the penetration depth distributions of events with 
and without shield hits has led to the conclusion that 91.3% ±3.7% of all rock 
events have shield tags. 

Our Monte Carlo simulation program tracks the decay products through the 
detector geometry and generates electronic hits in the same format as real 
data. The generator starts with a parent (or target) nucleon within a nucleus. 
The nucleon is considered to have a Fermi momentum chosen from the Bodek 
and Ritchie parameterization [[3~5j| . The spatial location and the atomic num- 
ber of the parent nucleus is chosen according to the composition and mass 
distribution of the detector. The rescattering of pions within a parent nucleus 
is generated according to a phenomenological model |16|]. Parameters of the 
model have been set to reproduce pion production by low energy neutrinos in 
z^-Ne and i^-deuteron reactions |ff7|. Inelastic intranuclear rescattering for K + 



mesons is not expected to be significant due to the absence of low-lying K + N 
resonances and is not simulated. The detector response used in this Monte 
Carlo program was verified against calibration data from the ISIS facility us- 



ing 7r, e, fj, and p beams at a variety of angles and energies[18 . 



3.2 Search for the Decay p — > vK + ,K + — > [i^v^ 



To choose the selection criteria and determine efficiencies we generated a large 
sample of K + — > [i^v^ Monte Carlo events from proton decay. Features of 
the typical Monte Carlo event shown in Figure [I] include a highly ionizing K + 
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which (usually) comes to rest, and emits a 236 MeV/c p + . The p + has an 
average range of 42 cm in the Soudan 2 detector. At the end of its range, the 
dE/dx of the muon is rising. The /i + comes to rest and decays into an e + , 
which gives an average of 3 hits. 

The results of the analysis for the proton decay simulation, the neutrino sim- 
ulation, the data, and the shield-tagged rock background are given in Table [I]. 
The simulated events are generated in the entire mass of the Soudan 2 detec- 
tor. Efficiencies within the fiducial mass are calculated by dividing the fraction 
of events which pass a set of cuts by the ratio of the fiducial mass to the total 
mass (0.79). 

The Monte Carlo events were first subjected to a simulated trigger, which 81 
events failed. Both Monte Carlo and data events were processed through a fil- 
ter program which applied containment criteria, reducing cosmic ray muons by 
a factor of more than 10 3 . 114 MC events, essentially all with hits outside the 
containment volume, were rejected. The remaining events were then scanned 
by physicists to remove the remaining non-contained events, mostly events 
starting or ending on module boundaries. At this stage, 204 p — > uK + ,K + — > /i + z^ 
events remained, compared with 367 data events without shield hits, 1008 
shield-tagged rock events, and 1923 events from the atmospheric neutrino 
Monte Carlo. The data corresponds to 3.56 fiducial kiloton-years of exposure, 
compared to 20.24 kiloton-years for the atmospheric neutrino MC. 

The required topology is two charged tracks with a common vertex. Events in 
which both tracks appear to be protons based on ionization and straightness 
are not included. These topology features were exhibited by 95 of the remain- 
ing 204 events from the MC. We also require that the K + candidate, which 
is usually the shorter track, has a length less than 50 cm. The muon track 
length distributions for all four event samples at this stage of the analysis are 
shown in Figure |[ The range of muons from the K decay is peaked at 43 cm 
whereas the background distributions are relatively flat. The muon range was 
therefore required to lie between 29 and 58 cm. Our final cut requires a visible 
muon decay electron having two or more hits. This requirement discriminates 
strongly against neutrino induced background, since the predominant back- 
ground is vn — > pTp and in our iron detector most pT are absorbed rather 
than decay after stopping. 

After all cuts, our efficiency for accepting p —> v K + ,K + — > p + v IJL events 
in the fiducial volume is 14%. Two atmospheric v Monte Carlo events pass 
all of the cuts and represent an expected background of 0.21 events, taking 
into account that we found the atmospheric neutrino flavor ratio {v^/v e ) to 
be only 0.61 of the expected value jn|. The rock event background in the zero 



shield hit sample is calculated to be 0.19 events. This was found based on the 
penetration depth analysis mentioned above, which found background in the 
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track data equal to 9.5% of the shield-tagged track sample. 
As shown in Table [3], one event in the data survives our cuts. 



MC Soudan2 Data 
Run 15878 Event 1786 
10-Dec-1989 16:25:57.93 




Fig. 1. Monte Carlo event of p — > v K + . The three projected views are shown. The 
xz/yz views correspond to the anode/time and cathode/time views. The xy view is 
based on matching anode and cathode hits using their time and pulse shape. For 
each hit, the pulse area is proportional to the recorded energy loss. The K + is the 
short heavily ionizing track on the left /left /right of the xz/yz/xy plot. Three hits 
from the e + decay of the /x + appear at the right/right/top. Scales are in cm. 



3.3 K + 



IT ' IT 







analysis 



We searched for the mode p — > v K + ,K + — > 7t + tt° by selecting events 
with a short heavily ionizing track (K + candidate), one other track and two 
showers. These features are illustrated in the Monte Carlo event shown in 
Figure 3. The short highly ionizing track is the K + before it decays. A tt is 
reconstructed from the two showers and the K + from the tt° and the second 
track, assumed to be a tt + . At these low energies the two 7's from the tt° are 
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Fig. 2. Muon range distributions (before the muon decay cuts are imposed) for a) 
236 MeV/c /i's from the p — > v K + Monte Carlo simulation, b) the data, c) the 
atmospheric neutrino Monte Carlo and d) the shield-tagged rock background. The 
numbers of events (including overflows) are the same as in the row labeled ll K range 
requirement" of Table 1. The shaded events pass the muon decay cut. 

usually identifiable. The K + mass is required to be in the range 100 MeV/c 2 < 
rriK+ < 660 MeV/c 2 . The ir + and 7r° momenta are required to be in the 
ranges 80 MeV/c< P?T + < 400 MeV/c and 40 MeV/c< p n o < 390 MeV/c. 
The invariant mass of the two shower system is required to be in the range 
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Table 1 

Numbers of MC and data candidate events for p — > v K + ,K + — > which 
survive the triggering, containment, topology, and kinematic cuts of this analysis. 
Events are generated in the full detector, while efficiencies in Table 3 are quoted for 
the fiducial volume. 

10 MeV/c 2 < m 27 < 290 MeV/c 2 . The resulting detection efficiency for this 
mode is 26%. 

The effects of these cuts on the proton decay simulation are shown in Table 2. 
No data events pass these cuts. The background from atmospheric neutrino 
interactions is estimated to be 1.05 events based upon our MC simulated v 
events. The background is from v e charged current and neutral current 
interactions, so the suppression factor of 0.62 is not used here. The estimated 
background contribution from rock events is 0.09 events, based on the one 
rock event which passed our cuts. 



3.4 Limit calculation 



Using the two K branching modes studied here we set a proton lifetime limit 
using the formula: 

r , _ N p xTx [e x x B X {K) + e 2 x B 2 {K)} 

-{p^uK + )> (1) 

B /il + jX 2 



Here N p = 2.87 xlO 32 is the number of protons in a kiloton of the Soudan 2 
detector, T = 3.56 kiloton-years is the fiducial exposure x B^K) are the 
detection efficiencies quoted above times the appropriate K + decay branching 
fractions in Table |3|. The /ij are the constrained 90% CL upper limits on the 
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Fig. 3. Monte Carlo event of p -> ^ K+; -> vt+vt . The if+ track of 

four heavily ionizing hits, ranges and decays at rest, yielding a pion track and two 
showers from the K + endpoint. The showers are the 7T° remnants; they appear in 
directions opposite to the pion track and are overlapping in the xz view. The xy 
view also shows the results of our track and shower fits. Scales are in cm. 

numbers of observed events, and are found by solving the equation 



0.10 



jni=0 Z^ri2=0 



P(n u 6i + /ii)P(n 2 , b 2 + /i 2 ) 



(2) 



with the constraint 



£l X B 1 (K) € 2 X B 2 (K) £•=! Si x 



^2 



E<=1 ^ 



(3) 
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Table 2 

Numbers of MC and data candidate events for p — > v K + ,K + — > 7r + 7r° which 
survive the triggering, containment, topology, and kinematic cuts of this analysis. 

where P(n, n) is the Poisson function, e _/i /i n /n!, and the 6j are the estimated 
backgrounds. With the one candidate observed for K + —> fJ^u^ and none 
found for K + — > tt + tt , the values of /i in Table |3] are obtained, and the 
combined limit without background subtraction is 4.3 x 10 31 years at 90% CL. 
With background subtraction, a limit 4.6 x 10 31 years is obtained. 
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Table 3 

Summary of detection efficiencies, branching fractions, backgrounds, and combined 
upper limits at 90% CL with and without background subtraction for the p — > v K+ 
channels analyzed. 
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4 Discussion 



The Kamiokande collaboration || has reported a limit for p — > v K + of 
10.0 x 10 31 years with 9 candidate events and an estimated background of 7.3 
events. Their unsubtracted limit is 5.1 x 10 31 years. 1MB has reported a limit 
without background subtraction of 1.0 x 10 31 years with 6 candidate events 
and an estimated background of 4.7. More recently, they report 5.0 x 10 31 



years without a background subtraction based on 14 candidates p0[. With 
background subtraction of 21.4 events, they calculated their limit to be 15.1 x 
10 31 years. Frejus has reported a background subtracted limit of 1.5 x 10 31 years 
with 1 candidate and an estimated background of 1.8. In all experiments, the 
largest contribution to the estimated background comes from quasielastic 
scattering. The existence of the atmospheric neutrino deficit then calls into 
question the reliability of background estimates in all of these experiments. It is 
not clear that rescaling the total rate to the observed rate leads to a correct 
estimate of the background for proton decay in this channel, since whatever 
is causing the atmospheric deficit may be energy dependent or different 
for and v^. We note the possibility of such differences by comparing plots 
2b and 2c. We would assign a resulting systematic error of as much as 100% 
to the background estimates for this channel. In view of this large systematic 
error we distrust the background subtraction in this channel and prefer to rely 
on the unsubtracted limit, which is unaffected by the background uncertainly. 
We note that the background level in this analysis, and thus the importance 
of background subtraction, is an order of magnitude lower than in the water 
cherenkov experiments. 

In summary, we have searched for the proton decay mode p — > v K + using 
two decay modes of the K + , K + — > and K + — > 7r + 7r°. We observe 

one candidate event for K + — > {i^v^ and zero candidates for K + — > 7r + 7r°; 
the estimated backgrounds are 0.40 events and 1.14 events respectively. Our 
combined lower lifetime limit at 90% CL is 4.3 x 10 31 years. Our limit with 
background subtraction is 4.6 x 10 31 years. 
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